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a^ 

^\ • The large deviation function obtained recently by Derrida and Lebowitz for the totally asym- 

metric exclusion process is generalized to the partially asymmetric case in the scaling limit. The 
O ■ asymmetry parameter rescales the scaling variable in a simple way. The finite-size corrections to the 

D ' universal scaling function and the universal cumulant ratio are also obtained to the leading order. 

00 

o 
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I . The asymmetric simple exclusion process (ASEP) is the simplest driven diffusive system where particles on a 

j^ I one-dimensional lattice hop with asymmetric rates under excluded volume constraints. Due to its simple but non- 
"t^ 1 trivial out-of-equilibrium properties, it has attracted much attention recently. We refer to W for a review of recent 
developments. 

For the prototype case of single-species, sequential updating dynamics, the time evolution operator of the probability 

distribution of particle configurations turns out to be the asymmetric XXZ chain PS| . The latter admits the Bethe 

rj^ ' ansatz solution for its eigenfunctions and eigenvalues when it is on a periodic ring. Due to its integrability, one can 

r^ I obtain many exact results of physical interest. In particular, the large deviation function (LDF) which describes the 

O ' distribution of the total current has been obtained recently for a ring of N sites with P particles under a periodic 

I ^i ' boundary condition ^g. The LDF also describes the height distribution of the Kardar-Parisi-Zhang (KPZ)-type 

growth models and is believed to be universal. To confirm the universality of LDF, Derrida and Appert B compared 

Vj ■ a cumulant ratio obtained from the analytic LDF with numerical simulations of several stochastic models believed to 

,__! , belong to the KPZ universality class. 

(^ ' Since the LDF has been obtained in |Q,|| for the totally asymmetric exclusion process (TASEP) where particle 

^^ hopping occurs only to the right, it would be desirable to calculate it for the partially asymmetric exclusion process 
^_^ [ (PASEP) where the particles can hop both to the right and to the left but with different rates. In this paper, we report 
Qv • on this generalization using the crossover scaling functions of the XXZ chain obtained previously in Q . Our method 
^\ , assumes from the outset that N is sufficiently large, but allows systematic evaluation of the finite-size corrections. 
We reproduce the universal scaling function of [^|j^ for the PASEP and find that the asymmetry parameter rescales 
the scaling variable in a simple way. We also evaluate the leading order finite-size corrections to the universal scaling 
G ' function and the cumulant ratio, 
i-rt This paper is organized as follows. In Sec. ||, we introduce the model and notation. In Sec. II], we make the 



G I connection between the present problem and the results of [^ and derive the LDF for the PASEP. The finite-size 
O ' corrections are evaluated in Sec. IV. Sec. [V| contains the summary and discussions, while Appendix shows the 



equivalence of two representations of the crossover scaling functions. 

II. MODEL AND THE LARGE DEVIATION FUNCTION 

We consider the dynamics of the one-dimensional model in a periodic lattice (ring) of N sites with P particles [E| . 
Each site j {1 < j < N) is either occupied by a particle {(Tj — —1) or vacant (ctj = -f-l). The PASEP considered 
in this work is defined by the following random sequential updating rule: During each time interval dt, each particle 
can hop to its right or left with probability ^(1 + e)di and ^(1 — e)dt, respectively, provided the target site is empty. 
e=l corresponds to the TASEP considered in |Q,g| and we work in the region < e < 1. Interpreting aj=±l as the 
local slope of an interface in (1 + 1) dimensions, one can map the model to the single step model Pfl], an archetype 
of the KPZ-class models. The quantity of main interest in this work is the total displacement Yt which is the total 



number of hops of all particles to the right minus that to the left between time and time t. In the single step model 
language, Yt is the total number of particles deposited between time and time t. 

Let a denote a system configuration {cti, . . . , a^} and Pticr) the probability of finding the system in a configuration 
a at time t. The master equation for the time evolution of Pt{(^) can then be written as 

^-^^^Y.^a\H\a')P,{a% (1) 

(T' 

where {a\H\(7') is the representation, on the basis where a^^ are diagonal, of the time evolution operator H given by 

H = -E{(^)<-.>1 + ^)-J-U + JK-I+1 - !)}• (2) 

Here, aj^ and a^ are the Pauli spin operators and (Tj=±1 are the eigenvalues of cr|. 

Next, following Mjm, we introduce Pt{a,Y), the joint probability that the system is in a configuration a and Yt, 
the total displacement, takes the value Y at time t, and let 

oc 

Ft{a;a)= ^ e"^P,(a,y). (3) 

F=-oo 

Then, Ft{(j;a) evolves according to 

'J:^^Y.(-\',\a')Ftia';a), (4) 

a' 

where 



^^ - E {^" (^) ->.>! + ^- (^) -7-U + iK-l+i - 1)} 



(5) 



The "Hamiltonian" —M is the asymmetric XXZ chain Hamiltonian studied, e.g., in g. Let \{a) denote the largest 
eigenvalue of M, regarded as a function of a. Then, one can show that 

(e"^*)=^i;^,(a;a)~e^(")* (6) 

as t —f oo and the long time behaviors of all cumulants of Yt are derived from A(q;). 

The LDF / describes the long time behavior of the distribution of Yt/t and is defined by 



f{y) = lim - In <( Prob 






(7) 



where v — limj_^oo(yi)/i is the mean current for a ring of finite size N . Note that v = dX{a) / da\ ^^q. This can be 
easily obtained from a first order perturbation calculation as 

N 
v^ep{l-p)N^;^-^. (8) 

Our definition of f{y) is slightly different from that of |jj^ in that we use the exact value of w, Eq. (g), in Eq. (J^) 
while @,| use its bulk value ep(l - p)N . Since (e"^*) - e^(°)* on the one hand, and {e'^^') = Ey=-oo Pi"ob[yt = 
F]e"^ ~ maxj, e*(-''(^)+"''+"'') on the other, the LDF is related to A(a) — av by the Legendre transformation 

fiy) = [^(") - "«] - ay: (9) 

y ^ ^[\{a) ^ av]. (10) 

Therefore, the largest eigenvalue \{a) of the asymmetric XXZ chain M determines the LDF. 



III. a(q) in the scaling limit 

In [HJq], A(a) for the case of e=l is obtained for arbitrary N and P. Then, one takes the scahng hmit, N — > cx), 



a — > 0, with the scahng variable aN'^^'^ and the density p = P/N fixed. In this scahng hmit, A(q;) takes the parametric 
form 



A(a) = aNp{l - p) + ^j P^-^f,/,(C) {e = 1), (11) 



where fk{C) are defined as 



a^2np{l-p)N^^h/^{C), (12) 



/3/2(c)=-5:L_^, (13) 

/5/2(C) = -E^' (14) 

n=l 



for \C\ < 1. To probe the region a^y2^Tp{l — p)N'^ < /3/2(~l), Ik{C) are analytically continued as 

fmiC) = -4Vi[- ln(-C)]i/2 - Y. i-^, (15) 

n— 1 

hMC) = ^^/^[- H-Cr^ - E ^' (16) 



for — 1 < C < 0, while for a-\/27rp(l — p)iV3 > /3/2(l), one may use the integral forms 

/•CXJ 

/5/2(C) = 27r / dsV^ln(l + Ce-"^). (18) 

Jo 

To generalize Eqs. (|ll| ) and ( |l2|) to the case of PASEP (0 < e < 1), wc limit our attention only to the scaling limit 
and use the results of Kim 0]. In [bI, the low- lying eigenvalues of the asymmetric XXZ chain near the stochastic 
line (a = in Eq. (|5|)) have been expressed as perturbative expansions in iV~^" with a scaling variable, which is 
essentially the same as aiV~^/^, held constant. Therefore the results of y] applied to the ground state energy (denoted 
as E^ in Pi) can be used immediately to obtain the LDF. 

The notations q, A, s, H, and i/ used in [^ translate into the present ones as p, (cosh a + esinha)^^, (sinha + 
ecoshQ;)/(cosha + esinha), {a + tanh^ £)/2, and tanh~ e, respectively. Using these and taking care of different 
normalization (A(a) = —E^j/A), one can rewrite Eq. (58a) of [g| as 

b„i,k ..n, . /7r\"/2 

n - 

m—1 k—1 

and Eq. (54a) of § as 



A(^)-^EE n-;V ^-(-H^) ' (19) 



iEEwH|!,,(,(^)--. 



m—l k—1 



In the above sums, Y^ for even m vanishes for the ground state and only odd-m terms are needed. For m odd, Yj'^^{z) 
with real z are defined as 






2i 



2t Jq e'"^ - 1 

The coefficients Xc and &„i,fe are recursively determined order by order in iV~^' ^ from a set of equations, as explained 
in 111 and Xc = -p/(l - p) + 0(iV~^/^), bm,k = b^f. + 0(iV~^/^). 6™^^ is the coefficient of x™ in the series expansion 
of (X]m=i '^■mX™') where a^ — a^m + 0{N~'^/'^) and the first few values of a^ needed in this work are given by 



a 







2p 



2 (l + p) 



2 ■ 



o_ / 2p 1 + llp + p^ 
"^ " V (1 - P)^ 18p(l - p) ■ ('') 

The eigenvalue expression Eq. (O) is a power series expansion in N^^^^ with the scaling variable aN^^^ > and 
e > fixed. (If a > and finite, the asymmetric XXZ chain is in the critical phase and hence the ground state energy 
and the low lying excitations possess finite-size corrections analytic in N^^.) When e ^ with another crossover 
scaling variable e^/N fixed, the infinite series Eq. ( [l9| ) reduces to a series in l/{e\/N). 

Inserting the zeroth order values of Xc and 6„i,fc, and keeping only the leading order terms in Eqs. ( [l9| ) and (|2q), 
one then obtains 



27riV3 f/ 47r2 



P(l-P) 



[A(a) - H = e I [-^-Y3(^)) ~ i^^Yii^))} (0 < e < 1), (23) 



ay/2T:p{l - p)N^ ^ (27rY]°(z)) . (24) 



Here v is the exact average current, ep(l — p)N^/{N — 1). The second term on the right-hand side of Eq. ( p3[ ) appears 
due to our choice of the exact v on the left-hand side of Eq. (p3|). Except for that, the similarity of Eq. ( p3| ) to Eq. ( pl| ) 
is obvious. One simply needs to relate Yj'^{z) to fk{C). In Appendix A, we show, by changing the integration contours 
of Eq. (|2^), that 2ttY^{z) is indeed nothing but a different form of /3/2(C'), provided the variables z and C are related 
by C = e'^^. So is —4tt^Y^{z)/3 of f 5/2(0) . Moreover, we show in Appendix that the analytic continuation of Eq. ( pl| ) 
to the region Im z > 1 naturally reproduces the analytically continued forms of Eqs. (|l5|) and ( |l6|) . Therefore, the 
generalization of Eq. ( pT| ) to e 7^ 1 is achieved by a factor e multiplying the right-hand side of Eq. (0) . Consequently, 
by Eqs. (|), (|o|), (ifl), and (M), one obtains the LDF in the form 



where the universal scaling function H{x) is given in the parametric form satisfying the relation 

„, s f 5/2(0) f 3/2 (C)- f 5/2 {C)f 3/2(0) 

nix) = ^=—, , (2b) 

V2f',/,(0) 

^ ^ /5/2 (O) - /3/2 (O) 

h/2(o) 

with ' denoting the derivative with respect to O. H(x), as defined here, is H(x + 1) of HJH], the difference originating 
from using exact v in Eqs. (0) and (|2^). Thus it has the following asymptotic behaviors 



x' 






for 


N 


I<1, 


2 
5 




^5/2 


for 


X 


-> oo, 


4 
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V^ 


:e|3/2 


for 


X 


— > — CX) 



H(x) ~ <^ _| ^/l 2;5/2 for:E->oo, (28) 

Eq. ( psi ) is the generalization of the result of 



IV. FINITE-SIZE CORRECTIONS IN DISCRETE DYNAMICS 

Finite-size correction is useful in comparing theoretical predictions with simulation data. In simulations, particle 
configurations are updated in discrete time steps, and to describe such situations, Eqs. (|l|) and ()^ should be replaced 
by their discrete time versions. For example, Eq. (0) is replaced by 

Fr+iia; a) - Fr{a; a) ^ ^ ^(a|Af |aV,(^'; "), (29) 

a' 

where one update interval is set as dt — 1/N and t — t/N. These difference equations reduce to the continuous time 
versions, Eq. Imj and Eq. (0), in the limit N ~* oo. Thus the leading terms in N of all quantities are the same in both 
versions. However, there appear differences in the finite-size corrections and we work in the discrete version. Using 
Eq. (H), Eq. (|) is then modified to {e"^') = Ea^t(f^;") ~ e'^^"''*' where 

M«) = iVln(^l + ^). (30) 

Therefore the LDF is the Legendre transformation of fi{a) — av. 

From Eqs. (|l9|), (EOT), and (pG), fJ^ia) — av is written as, including its next leading term. 



2^^' -[^(a) - av] ^ 6{/5/2(C) - fy,{C)} - ^J P^^^f„^{Cf . (31) 



The last term on the right-hand side of Eq. ( pl| ) arises from the first nonlinear term in the expansion ^{a) — 
\{a) — \{a)'^ /2N + ■ ■ ■. The leading correction term in iJ,{a) — av is of order iV~^/^, while that in \{a) — av is of 
order N^^. Since the leading correction to ay^2TTp{l — p)N^ is also of order N^^, using Eq. (p3) and Eq. (Bl|), one 
finds that 



ny)-Jp^'~'^ 



^ep{l-p)) V TTiV Vep(l-P)' 



(32) 



where Hi{x) is determined from Eq. ( p7| ) and 

rj , ^ /3/2(C) ,„„x 

Hi{x) = '— . (33) 

The correction term shows dependence on the particle density and the asymmetry parameter, and hence is not 
universal. The asymptotic behaviors of Hi (x) are 

-2x^ for |a;| < 1, 

Hi{x) ~ { -3x^/{2n) for x ^ oo, (34) 

-27r|x| for x — > — oo. 

Another quantity of interest concerning the finite-size correction is the cumulant ratio considered in [p| . It is defined 
as 

/y3\2 

lim i?, = lim -4-^-04-, (35) 



where (Y^)c are the cumulants of Yt and are evaluated from 



Urn 

t — *oo t 



(y,"), _ dXa) 



da"^' 



(36) 



a=0 



Using Eqs. (EJ) and (M), we find 



t— >oo t Z 



l_2e,/^ii^ + 0(7V-i) 



(y*'>c 



lim 

t-^oo t 



eN^W-p)?Al-^-f 



lim 



{Yt 



^=e7V«/2[p(l-p)]5/V/2(l^ + 



t^oo t irv r/j 12 2 

Therefore the cumulant ratio has an 0{N^^i'^) correction term; 




15 9V2 



8V3 [l + 0{N-^)\ 



(37) 



lim Rt 

t — 'OO 



3 _ sVa V 

2 <3 ) 



^ + 9^ _ 8V3 



l + 2.J^ii^ + 0(^-) 



(38) 



We note in passing that in the continuous time version, our method shows 



{y? 



lim 



eN^/M^-p)f/^^ 






(39) 



This is in exact agreement with the expansion derived, with the help of Stirling's formula, from Eq. (6) of Derrida 
and Mallick . 



V. DISCUSSIONS 



The main results of this paper are Eqs. (|25D, (32), and (|38|). The universal scaling function of the LDF, H{x), first 
defined in HM| for the TASEP, is reproduced for the PASEP in Eq. (p5|). The only change in this generalization is 
the modification of the scaling variable by a simple factor e, the asymmetry parameter. Physically, this is equivalent 
to a rescaling of time by e. Non-trivial e-dependence of A(a) appears only in higher orders of iV~^'^ in Eq. (M). 
To compare analytic results with simulation data, the finite-size correction terms in the discrete time dynamics are 
important. They are derived for the LDF and the cumulant ratio in Eq. ( p^ ) and (p8|), respectively. One sees that 
the finite-size corrections in the discrete time dynamics are of 0{N~^''^). Also they depend on p and e explicitly in 
both versions implying that they are non-universal. 

Instead of the statistics of Y^, the total displacement, one could have asked for the statistics of the displacement 
across one bond. In this case, one has to deal with an asymmetric XXZ chain with a twisted boundary condition. 



' N+l 



?^°'at and a 



N+l 



af, and analysis similar to that presented here can be carried out Iq]. In particular, if 
Jf is the displacement across the A^-th bond, one can show that (e""'*) ~ (e"^'/^) ~ eH^'/^)'^ (in the continuous time 
notation). Therefore, the LDF and the long time behaviors of the cumulants of Jt are the same as those of Yt/N. 
This is why Eq. (BOh agrees with the result of JTl where liTat_,^{Jf)c/t is obtained. However, {Jt^)c — ((^t/A^)^)c, the 
surface width in the growth model language, saturates to a finite value of 0{N) as i — > oo. 
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(b) 

FIG. 1. Contours for 7i (a) and I2 (b). Re z is assumed to be negative and [Im z\ < 1. The semicircles about the origin 
have the small radii 5, which are set to be zero in the last step. 

APPENDIX A: PROPERTIES OF Y^{z) 



In this Appendix, we show the equivalence of fi+rn/2{C) and Y^{z) {m = 1, 3, 5, . . .). The former will be defined 



later extending the definitions of /3/2(C') and f5/2{C), and the latter is defined in Eq. (21). We take Eq. ( |21| ) as 
defining Y^{z) for any complex z. 



1. Simple form of Yl\z) 



We first pay attention to Y-[^(z) since Y^{z) (m — 3,5, . . .) can be evaluated from Y-[\z) through the recursion 
relation, dY^^2{z)/dz = -{m + 2)Y^{z)/2 |]. Eq. ^ for m = I is written as 



0..O. ^ V .M/2 2, .,3/2 r Az + i-tf/^-iz + i + t)^/-' 
2Y^{z) = -i{z + iY'^ + -{z + if/^ + I dt^ ' ^ ' 

•J Jo 



i/2_r;,_i , ^u/2 



^■( -^1/2^2 .,3/2 , r . M-'-^r-i^-' + i) 

+i{z-i)' +-{z-i)' +1 dt e'^t -1 



(Al) 



Suppose |Im z| < 1 and let /i and I2 be the first and the second integrals in Eq. (Kfl), respectively. The two 
integrations are over the positive real axis of the complex-t plane, denoted by K in Fig. 111. 

Our method is to deform the contours in the complex-f plane such that only simple integrals remain and the additive 
terms cancel out. Each integral has two terms. For the first term of /i, K is deformed to Ai + A2 + A^ as shown in 
Fig. |l|(a) while for the second term of /i, to Bi+ B2+ -B3 in Fig. |l|(a). Similarly, K is deformed to Ei+ E2 + E3 and 
Di + D2 + D^ for the first and second terms of /2, respectively, as shown in Fig. p](b). We then have 



h 



dt 



A1+A2+A3 



jz + i-t) 



1/2 



dt 



B1+B2+B3 



(z + i + t)^/^ 






., ,«.l(i±i)!^-r"-"'*l(i±i)!^ 



ds 



g7r(z+i){ _ 2 



gl/2 



i{z + iY'^--{z + if/^+ / ds , ^., 



g-7r(z+i)g7rs _ ^ 
(-S)V2 



ds- 



,1/2 



, — 7r(2+'i)g7rs ]^ ' 



(A2) 



-i 



Di 
-z + i 



(a) 



(b) 



FIG. 2. Contours for /i (a) and I2 (b). Re z is assumed to be negative and 1 < Im z < 3. Compared with Fig. hi, a pole is 
placed inside the contour for /i and the integration over the semicircle about the origin makes a different (sign-changed) value 
in h- 



where 9 — Arg(z + i). Similarly, 



/2 = / dt 



(z-i-t) 



1/2 



<-\,^i^^^ 



l^^(,_,)3/2(l_^)l/2 



d9 



dt 

D1+D2+D3 



jz-i + t) 



1/2 



ds 



g7r(z— 0? — 1 
g7r(z— i)g7rs _ ^ 



10 

ds 



1 _(^^^)3^2(^_^ 
^ g-7r(2-i)4 _ I 

5I/2 



-^{z-^y/'-^{z-^f/' + 



g-7r(z-i)girs _ J^ 
(-.)V2 



ds 



T{z-i), 



1 



ds- 



,1/2 



, — 7r(2: — z)g7rs ^ 1 ' 



(A3) 



with 6* = Arg(-2;-|-i). 

We note that the branch cuts for the square-root functions in /i and I2 are in the opposite directions H, so the 
two integrals having the factor (— s)^'^ in their integrands cancel out when /i and I2 are added. Therefore we arrive 
at the conclusion that 

/•oo 1/2 

^iW=/ ds _^^ ^ (|Imz|<l). (A4) 

Next, consider the region 1 < |Im z| < 3. If 1 < Im 2: < 3, the contours shown in Fig. |^ change to those shown in 
Fig. g. Compared with Fig. |l], a pole at t — 2i is placed inside the contour for /i and the direction of the integration 
over the semicircle about the origin is reversed for /2. These changes produce extra contributions 2i(z — i)^" for both 
/i and /2- Therefore, we obtain 



r°(z) = 2i(z-i)^^^ + 

= -2(-z+i)i/2 + 
Similarly, for — 3 < Im z < — 1, we find 

Yi°(z) = -2i(z + i)i/2-f 

= _2(-z-z)i/2 + 



ds- 



«l/2 



ds- 



,1/2 



(1 < Im z < 3). 



(A5) 



ds- 



U/2 



ds- 



,1/2 



(-3 < Im z < -1) 



(A6) 



When |Im z\ increases further, more and more poles are placed inside the contours for Ji and /2, and corresponding 
residues should be added. But recalhng that Y^^{z) is used to express real a and A(a), Eqs. ( |A^ ) and (A5) are 
sufHcient for our purpose. 

2. Relation bet'ween f\+m/2{C) and Y^{z) 

We now make the identification C = e"^ with C real. C > if z is real, and — l<C<Oifz = ~x + i^ with a; > 0. 



In both cases, we have, from Eq. (A4), 



yO(z) = / ds—i- 

^^ ' Jo C-^e-^" + 1 



Eq. (A7) admits a series expansion 



_. oo 



(-l)"+iC 



27r ^ — ' n 

n=l 



3/2 



(A7) 



(A8) 



when |C| < 1. The second branch in the region — 1 < C < is obtained if z = —x + z+ with x > Q. In this case, using 
Eq. ®, 

9 r°° «i/2 

r°(z) = -4[-ln(-C)]V2+ / 

V^r Jo 



-:i-H-cr'+l± 



(-l)"+iC" 



27r ^^ n3/2 

n—1 



(A9) 



Comparing Eqs. (A7) and (A9) with Eqs. (Of) and mW, we have 

/3/2(C) = 27rY^\zl 

provided C = e'^^, the two branches of —1 < C < corresponding to |Im z| = 1^ and |Im z\ — l'^ , respectively. 
Next, we define fi+,n/2{C) {m = 1, 3, 5, . . .) by Eq. ( AlC ) and the recursion relation 



(AlO) 



d(lnC) 



/l+(m+2)/2(C) — /l+m/2(C), 



with the initial condition /i+to/2(0) = 0. Then /i+m/2(C') takes the form 



/l+m/2(C) = 2^ ^l+m/2 
n—1 



in the first branch and 



'-A^[-H-c)r'^ + Y, 



(-l)"+iC" 






(All) 



(A12) 



(A13) 



in the second branch. Comparison of the recursion relations of Y,^{z) and fi+m/2{C) then leads to the identification 

h+rn/2{C) = (-7r)(™-i)/227r|^r^(z). (A14) 

For example, fy2{C) = 27rYi"(z), f^^^iC) = -A7r^Y^°{z)/3, h/2{C) = 87r3y5"(z)/15, etc. 
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